Gapped Two-body Hamiltonian whose Unique Ground State is 
Universal for One-way Quantum Computation 
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Many-body entangled quantum states studied in condensed matter physics can be primary re- 
sources for quantum information, allowing any quantum computation to be realized using measure- 
ments alone, on the state. Such a universal state would be remarkably valuable, if only it were 
thermodynamically stable and experimentally accessible, by virtue of being the unique ground state 
of a physically reasonable Hamiltonian made of two-body, nearest neighbor interactions. We intro- 
duce such a state, composed of six-state particles on a hexagonal lattice, and describe a general 
method for analyzing its properties based on its projected entangled pair state representation. 

PACS numbers: 03.67.Lx, 03.67.Pp, 42.50.Ex 



Many-body entanglement is fundamental to the under- 
standing of complex condensed matter systems, as well 
as a primary resource for quantum computation. A sur- 
prising result in quantum computation is that certain 
entangled states can be employed to perform arbitrary 
quantum information processing tasks, merely by sys- 
tematically measuring single sites in different bases, in 
a method known as "one-way" quantum computation [![. 
If such universal resource states are available, this ap- 
proach potentially simplifies experimental requirements 
by employing only measurements, and not multi-qubit 
gates normally needed. Exotic physical properties may 
arise in these states due to their large amount of entangle- 
ment; many methods have been developed in condensed 
matter theory to study such systems, including the ma- 
trix product state formalism Q or more generally the 
projected entangled pair (PEPS) representation of states 
1- 

The special entangled states which make arbitrary one- 
way quantum computation possible unfortunately seem 
to be difficult to realize naturally. Ideally, such uni- 
versal resource states could be obtained as the unique 
ground state of a naturally occurring Hamiltonian, one 
with only nearest-neighbor two-body interactions. If this 
were the case, especially if an energy gap existed between 
the ground and first excited states, the one-way quantum 
computation could be robust against quantum noise and 
decoherence of the entanglement. However, no such par- 
ent Hamiltonian exists for any of the presently known 
resource states of one-way quantum computation. 

Many efforts have been made to construct the desired 
many-body entangled state such that it could be the 
ground state of a naturally occurring Hamiltonian. The 
first and best known resource state is the cluster state, a 
simple entangled state on a two-dimensional square lat- 
tice; unfortunately, it cannot be the exact ground state 
of any naturally occurring Hamiltonian Perturbative 
approaches providing a Hamiltonian whose ground state 
approximates that desired have been developed 
A nice scheme for constructing universal resource states 



has been proposed and has yielded many interesting 
examples [8]. Based on this, a mixed approach can be 
taken, using a ID Hamiltonian to create chains, that 
are then coupled by two-body unitary operations 0, §] 
to form a 2D resource state. Matrix product state 2j 
techniques allow any measurement of these ID chains to 
be computed efficiently, on a classical computer, how- 
ever, implying that they alone are insufficient for quan- 
tum computation. Two-dimensional many-body entan- 
gled states are thus likely to be essential for arbitrary 
quantum computations, but few techniques are presently 
known for finding local 2D Hamiltonians with the desired 
ground states. Properties of such states generally remain 
intrinsically hard to determine 

Here, we present results from a new approach to 
studying the quantum informational and physical prop- 
erties of 2D many-body entangled states using the PEPS 
representation. On the one hand, this representation 
naturally includes many-body entanglement in its state 
description Q and hence facilitates understanding of one- 
way quantum computation schemes[il.[lol|. On the other 
hand, methods have been developed to study the physi- 
cal properties of PEPS states as ground state of parent 
Hamiltonians |ll|. Combining these insights, we are able 
to construct the first example of a system which is both 
the unique ground state of a gapped two-body nearest- 
neighbor Hamiltonian and a universal resource state for 
one-way quantum computation. 

Building on PEPS - A good place to start in con- 
structing the desired state is with an example which il- 
lustrates the challenge, based on the well-known cluster 
state. Consider the state IV'peps) defined on a square 
lattice (Fig.[T^) where each pair of nearest-neighbor sites 
are connected by singlets \<p) = |00) + |01) + 1 10) - 1 11) 
(suppressing normalization). On sufficiently large lat- 
tices, starting with IV'peps)' anv quantum circuit can be 
efficiently simulated by measuring all four qubits at each 
site (on the boundary, two or three) in appropriate time 
sequences and measurement bases [lfj. 

IV'peps) 1S * ne unique ground state of a gapped two- 
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body Hamiltonian, as it is simply a tensor product of 
two-body entangled states. However, the multi-particle 
measurement required to make this state universal flpj| 
is generally disallowed in one-way quantum computation 
models. Still, if the four qubits at each site were treated 
as a single 16 dimensional particle, the model could be 
interpreted as giving the desired result, a universal re- 
source state for one-way quantum computation with a 
gapped two-body parent Hamiltonian. And while use of 
16 dimensional particles is experimentally unrealistic, the 
idea of using a description in terms of singlet pairs does 
provide a good starting point for constructing simpler 
states. 
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FIG. 1: Projected Entangled Pair State (PEPS) represen- 
tation of 2D (a) square I^peps) an d (b) hexagonal I^peps) 
lattice states. Filled circles connected by solid lines denote 
virtual singlet pairs \ip). Dashed circles denote projection of 
virtual qubits into physical states; red and green ones corre- 
spond to sublattices A and B in hexagonal lattice respectively. 



Specifically, consider the set of states given by pro- 
jecting lattices of singlets into smaller subspaces. For 
example, the projector Poster = 1 6) (0000 1 + |1)(1111|, 
applied to all sites of the square lattice state gives 
the cluster state on a square lattice [l(J, {^cluster} oc 

Pciusterl^PEPs) ' wnere 1 6) and |1) are the physical qubits 
in the cluster state model. In this "PEPS representation" 
picture, the physical PEPS state is defined by two ele- 
ments, a lattice of "virtual" singlets (connecting neigh- 
boring sites in the lattice), and a set of projectors which 
act on lattice sites. Not all PEPS states are universal for 
quantum computation; only a few, such as |\& cluster), are 
known to be universal. 

Compared with IV'peps) w ith 16-dimensional particles, 
[® Cluster) employs only qubits at each site, and hence is 
more experimentally accessible. Unfortunately it cannot 
occur as the exact ground state of nearest-neighbor inter- 
actions jj], and its parent Hamiltonian involves at least 
five-body interactions. Moreover, it is known that PEPS 
states generally disallow low dimensionality particles si- 
multaneously with short interaction ranges 11[. Never- 
theless, this line of thought, using PEPS states, can in- 
deed lead to a universal resource state which is the unique 
ground state of a gapped nearest-neighbor Hamiltonian, 
while also being composed of particles of relatively low 



dimension, as we now show. 

The tri-Cluster State - The structure of the lattice 
of singlets, and the choice of projectors, in the construc- 
tion of PEPS states, provide powerful degrees of freedom 
for exploring interesting new states. Two specific insights 
from the above examples illustrate this freedom: 

1. Instead of on a square lattice, a cluster state defined 
on a hexagonal lattice \^ cluster) IS a ^ so universal [12j. On 
a hexagonal lattice of singlet pairs (Fig.QJ)), the projector 
defining this cluster state is Pester = |0><000j H-| 1> <111| , 
giving |*g£f sier ) oc PgtZterlfpKPs)' where the labels de- 
note left-right-up and left-right-down virtual qubits on 
sites in sublattices A and B, respectively. 

2. An alternative projector can be chosen: P' = 
|6)(100| + 1 1> <01 1 1 or P" = |0)(010| + 1 1> <101 1 : these re- 
sult in PEPS states different from |*c^ ter ), but onl y 
by local Pauli operations. Hence, a modified local mea- 
surement scheme still exists, allowing these states to also 
be universal. 

We now introduce a new state, the tri-Cluster state 
I 'I' trie), which is motivated by these two insights, and 
has properties we desire. This is defined in the PEPS 
representation on a two-dimensional hexagonal lattice 
(Fig. [TJd), with projectors 



PtriC = |0)(000| 
+ |2}(100| 
+ |4)(010| 



|1)<1H| 
|3)(011| 
15X1011, 



(1) 



using the same labeling scheme as above, such that 
Y&tric) oc PtriC I^keps) ■ Hence, at each lattice site there 
lives a 6-dimensional particle. 

Intuitively, \^>tric) is universal because it is closely re- 
lated to the standard cluster state. Specifically, \$>tric) 
projected into the subspace spanned by {|0), |1)} is the 
same as cluster), as are also the states given by \^tnc) 
projected into {|2),|3)} and {|4),|5)}, up to local Pauli 
errors. Thus, \^tric) is like a "superposition" of three 
cluster states. Computational qubits are encoded in the 
virtual qubits and operated upon by measuring the phys- 
ical particles. Although the three subspaces of \^tric) 
cannot be decoupled physically, they may be employed 
independently in processing encoded qubits with a suit- 
able choice of measurement basis, as detailed later. 

The most interesting nontrivial feature of \^tric) is 
that it is the unique ground state of a gapped two-body 
Hamiltonian, and we begin with that. 

Uniqueness &; Gap - The fact that \^tric) occurs 
as the unique ground state of a gapped two-body Hamil- 
tonian is very surprising, as on the one hand the ground 
states of two-body Hamiltonians are rarely exactly known 
and on the other hand simply constructed states do not 
usually have simple parent Hamiltonians. Even the one- 
dimensional cluster state requires 3-body interactions 
in its parent Hamiltonian. Below, we give a two-body 
nearest- neighbor Hamiltonian H tr ic which has \^tnc) 
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as its ground state. Furthermore, we prove that \^>tric) 
is the only ground state of Htric and the Hamiltonian 
has a constant gap independent of system size. 

The central step in constructing H tri c and studying 
its properties is to find the support space S a b of the re- 
duced density matrix of any two nearest-neighbor par- 
ticles a and b in the state (a, b are in two sublattices 
A, B respectively). This is accomplished by first finding 
the corresponding support space S^ b EPS of the six virtual 
qubits on site a and b, in the PEPS picture, and then 
computing S a b <x PtricS p b EPS - For example, when a is 
to the left of b (Fig. [2j, virtual qubits 1 to 6 on those 
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FIG. 2: One representative site with particles a and b, and 
neighboring boundary, in the hexagonal lattice of IV'reps)- 
Filled circles connected by solid lines represent virtual singlets 
\<p) and dashed circles indicate sites projected to obtain the 
physical state. 

sites are only connected to virtual qubits a, (3, 7, S else- 
where. By tracing out a to S from the 5 singlet pairs, 
we find S^ b EPS for virtual qubits 1 to 6 to be spanned 
by |±) 1 |±) 3 |vj)m|±)6|±)6, where |±) = (|0) ± \l))/y/2 
and \ip) is the singlet state. This 16-dimensional space 
is then projected to give S a b for the depicted lattice site. 
S a b is different for the three bond directions in a hexago- 
nal lattice, i.e. a to the left of, to the right of, and below 
b. 

Providing a two-body Hamiltonian with \^tric) being 
a ground state is straightforward. The Hilbert space of 
two neighboring sites a, b is 36-dimensional, larger than 
the dimension of S a b- Therefore we may choose any non- 
negative Hermitian operator h a b on the two sites that has 
S a b as its null space, such that h a b\^ trie) = for every 
h a b- Thus, \^tric) is a ground state of the two-body 
Hamiltonian Htric — S & hab, where the summation is 
over all nearest-neighbor pairs. However, the key is to 
construct H tr ic such that \^tric) is the unique ground 
state, and it turns out the above procedure does work. 

Specifically, let h a b be the projection operator h p ab 
which projects onto the 36 — 16 = 20 dimensional sub- 
space orthogonal to S a b, giving the total Hamiltonian 



HtriC — 



E 



(h P 



b + Ka 



hi 



(2) 



The summation is over sites a in sublattice A and the 
three terms h? ab , h ba , h p b correspond, respectively, to three 

a 

bond directions where a is to the left, to the right and be- 
low b. The Hamiltonian is hence invariant under transla- 



tion along sublattice A. An explicit expression for H tr ic 
in terms of spin operators can be given (see appendix). 

The specific H tr ic we have presented has \^>tric) as 
its unique ground state. This is shown by verifying the 
condition [ll| that for any region R of spins in I'&tric), 
the support space Sr of the reduced density matrix on 
R satisfies 



Sr 



f) S ab 

(ab) 



l R\ab ■ 



(3) 



where the intersection is taken over all neighboring pairs 
ab and In\ a b is the full Hilbert space of all spins in re- 
gion R except a and b. For every possible configuration 
containing three or four connected sites in \^tric) the 
condition is confirmed by direct calculation. To check 
the condition for larger regions, it is useful to notice that 
any region in \^>tric) containing more than one site is 



injective [llj|. By Lemma 2 of [llj, 1. if regions R\ and 
i?3 are not connected and R2 and i?3 are injective, then 

Sr 1 ur 2 ur 3 = (Sr 1 ur 2 ® Ir 3 ) n {Sr 2 ur 3 ® iRi) 2- if re- 
gions R\, i?2, R3 are all injective, then Sn lU ]i 2U R 3 — 
(S Rl ur 2 ® Irh) H (Sr 2 ur 3 ®I Rl )n(S Rl uR3 ®Ir 2 ). Hence 
for a region R containing more than 4 sites in \^tric)i Sr, 
is the intersection of all four-body support spaces in R. 
By induction, it follows that condition Eq.([3]) is satisfied 
on Y&tric) for any R. Therefore, Y&tric) is the unique 
ground state of Htric- 

Htric is also gapped; an energy gap 77 above the ground 
state exists, which is constant as the system size goes to 
infinity. The existence of this gap guarantees protection 
of l^tric) against thermal noise, independent of system 
size. 77 can be bounded. First, we show that r\ is greater 
than A, the gap of another Hamiltonian K which also 
has Y$tric) as its unique ground state, but has four-body 
terms instead of only two-body terms. We then bound A 
above a positive constant value. 



FIG. 3: Regrouping of lattice sites in tri-Cluster State into 
disjoint blocks, each containing two sites. 

Consider the Hamiltonian K for a re-labeled ver- 
sion of \^tric), i n which particles are regrouped into 
disjoint blocks each containing two nearest-neighbors 
(Fig. [31). Let K = J2m n kmn, where m,n denote 
two connected blocks, each containing two particles 
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m^mM and n''',n' r ' respectively, and k mn is pro- 
jection onto the orthogonal space of the four-body 
reduced density matrix on , , , (assum- 



ing mM and are connected). 

J2ab^ab - iEtim (' l I[^- M +' lP 



Then i? tri c 



> 



' ^[!] n [r]) 



.6—4 L~*mn 

and fc mn are both non-negative operators with the same 
null space, so the last inequality holds for some posi- 
tive number [i. Assume that the gaps of the projec- 
tors h p ab and k mn are both 1. Direct calculation gives 
/j = |. As discussed in the uniqueness proof, K also 
has \^tric) as its unique ground state. Using this, we 
find i] > i/iA = |A. The gap A can be bounded 
by showing that K 2 > cK for some positive constant 
c K 2 = i^2 mn k mn ) 2 = K + ^2 mnm , n ,(k mn k m i n r + 

krn'n'k'mn) ^ K F '^^tnimrij k 7nni k nln ^ -j- k nln jk mni7 Ti 7 

and rij are blocks connected to m. The last inequal- 
ity holds because when region mn and region m'n' do 
not intersect k mn k m > n > + k m > n /k mn > 0. Direct calcu- 
lation shows that (Fig. k mni k mn] + k mnj k mnz > 
for = (1,2), (1,3), (2, 4) or (3,4) and 



> 



-±k 

3 ^mrii 



for = (1,4) or (2,3). 



Summing over all consecutive 



m, 



and 



gives 



.mnkmn- There- 



fore -ftT 2 > \K , giving A > \. Finally, we find a lower 
bound on the gap 77 of H tr ic of rj > g A > 

Universality - \^tric) is a universal resource state, 
because of properties it inherits from the cluster state. 
Similar to a cluster state, computational qubits are en- 
coded in the virtual qubits, and the active computa- 
tional state flows along the lattice as measurements on 
the physical states are performed. In contrast, however, 
with \^tric) extra Pauli errors occur, thus necessitat- 
ing additional analysis. Below, we describe the different 
steps necessary, focusing on initialization and readout, 
one-qubit gates, and a two-qubit gate sufficient for uni- 
versality. 

Initialization and readout: Just as with the cluster 
state, with \^tric)i measurement in the six-state basis, 
{|0) • • • |5)} accomplishes several tasks. First, such mea- 
surement detaches unnecessary sites from their neighbors 
(up to a known Pauli error). Next for state initialization, 
it gives a post-measurement state with an encoded qubit 
projected into |+) (when the outcome is 0, 3 or 4) and 
|— ) (for outcomes 1, 2 or 5). At the end of computation, 
the encoded qubit can also be read out in this way, giving 
(for 0, 2 or 5), and 1 (for 1, 3 or 4). 

One-qubit gates: Similar to gate implementations with 
the cluster state, once a line in the lattice has been de- 
tached from the rest, appropriately measuring a site in 
the line performs a single qubit rotation, up to a known 
Pauli error. Specifically, measuring in the basis {|0) ± 
e' ie |l), 1 2) ± e* e |3), |4) ± e" ie |5)} implements operation 
{HZ(0), XHZ{9),ZHZ{9),YHZ{0),ZHZ{0),YHZ(0)}, 
respectively, on the encoded qubit (using standard nota- 



tion for qubit gates, with Z(9) denoting a rotation about 
z by angle 0), up to pre-existing Pauli frame errors from 
detaching the line. 

Two-qubit controlled-Z gate: Measurement of two 
vertically connected particles a and b implements the 
final ingredient needed for universality, a controlled- 
Z gate CZ a i„, just as with the cluster state scheme, 
but with some additional Pauli frame errors. Specifi- 
cally, measuring in basis {0...5} = {|0) ± |1), |2) ± 
1 3), 1 4) ± 1 5)} implements the two-qubit operation 
(X^Z^H a ) (X?Z?H b ) ® (X^X^CZ ab X^X^) on 
the two adjacent encoded qubits. For x € {a, b}, u x = 1 
for x measurement outcomes 1, 3, or 5; v x — 1 for out- 
comes 2, 3, 4, or 5; w x = 1 for 4, 5; and u x ,v x ,w x are 
otherwise. Much like for the cluster state, when embed- 
ded in a larger circuit, more complicated configurations 
arise in implementing a controlled-Z gate (see appendix), 
but the principles of propagating a Pauli frame remain 
the same. 

Conclusion - Y&tric) is a remarkable entangled many- 
body state which is both universal for one-way quantum 
computation and the unique ground state of a gapped 
Hamiltonian H tr ic, made of local two-body terms. While 
imperfect, due to use of six-state spins, it steps far closer 
to physical realizability than previous models. Moreover, 
the methods introduced here, based on the PEPS repre- 
sentation, arc very general. These analysis methods lead 
directly to a number of additional universal states, and 
deepen connections between the study of complex con- 
densed matter systems and quantum information science. 
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Appendix 



B: Two-qubit gate operation: example 



A: Explicit Hamiltonian with spin operators 

Y&tric) is the ground state of a Hamiltonian Htric, 
which can be expressed explicitly as being a sum over 
two-body nearest-neighbor interactions between six-state 
spins on a hexagonal lattice. Below, we give such an ex- 
plicit expression with spin operators, as Hf riC . Let the 
spin operators for particles a and b along directions x, 
y and z be S aa , S ay , S az and S bx , S by , S bz respectively. 
S+ = S x + iS y and 5_ — S x — iS y are the raising and 
lowering spin operators. H£ riC is translationally invari- 
ant along sublattice A and is a sum over three sets of 
local terms at every site a in A: 



(hab 



l>ba 



(4) 



h a b, hba and h b describe interactions along the three bond 

a 

directions at each site a, where b is to the right, to the 
left and above a respectively. They can be expressed 
explicitly as: 

h a b = 

2(2S a , - 5)(2S Uz - 3)(2S az - 1)(25„, + 1)(4S„, + 11) 

(2S bz + 5)(2S b , + 3){2S bz - 1)(2S 6 , + 1)(4S 6 . - 11) 
- 7bV^S a+ (2S az - $){2S az + 3)(2S az - l)(2S a , + 1) 

(48S£ + 64S£ - 280S£ - 2725^ + 67) 
+ 75^2(48^ -MSI - 280SI +2725 0a +67) 

S b+ (2S bz -5)(2S bz - 3)(2S b "- l)(2S bz + 3) 
+ 4VWSl + (2S az - 1)(2S„, - 3)x 

(128S£ + 560S£ - 2840S£ - 3M8S bz + 675) 
+ 4VTo(128S| - 5605* + 2840^ - 38485a, - 675) 

Si (2S bz -5){2S bm -3,) + h.c 



The hexagonal lattice structure of \^tric) can com- 
plicate implementation of controlled-Z gates when em- 
bedded in larger circuits. For example, if computational 
qubits are encoded only in every other line, then extra 
steps arc needed to bring two qubits together for the two- 
qubit gate. Below, we illustrate how this can be done, 
analogous to the complicated case presented in the orig- 
inal cluster state paper [![. 

Consider the configuration shown in Fig. |4l and sup- 
pose information is made to propagate from left to right 
and the top and bottom lines connect only through sites 
c and d. The initial state |$ a fc) is input at sites a and b. 
Red sites e and / are detached from the rest by measur- 
ing them in the six-state basis, {|0), . . . , |5)}, giving mea- 
surement results j e and jf. The different measurement 
results obtained contribute different Pauli errors to the 
two-qubit operation done next. Measuring the green sites 
a, 6, canddin basis {6.. .5} = {|Q}±|1), |2)±|3), |4)±|5)}, 
gives results i a , i b , i c , id, and performs the two-qubit op- 
eration operation H a H b ®CZ ab on |$ a f>) up to some Pauli 
operation, leaving the final states output at sites g and 
I. 
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where h.c. denotes the Hermitian conjugate, as usual. h ba 
can be obtained by exchanging a, b in the above. h b is: 



-25{2S az - 5){2S az - 3){2S az + 3)(2S az + 5) 
25S a a+ (2S az -5)(2S az -l) 

(224S| - 16S£ - 1968S£ + 40S£ + 3550S bz - 9) 



(4165^ - 80S£ 
+ h.c. + (a <^ b) , 



3600S£ 



520S? 



599AS h 



125) 



where (a b) denotes an exchange of a and b in the pre- 
ceeding expression. Because each positive-semidefinite 
local term in H* riC has the same null space as h p ab , H£ r 



triC 



is gapped and has \^tric) a s its unique ground state for 
the same reasons as for H tr ic- 



FIG. 4: Measurement pattern for simulating control-Z gate. 
a,b are input sites and g, I are output sites. Green sites are 
measured in basis {|0> ± |1), |2) ± |3), [4) ± |5}}. Red sites e, / 
are measured in computational basis {|0), . . . , j5)} and their 
results may also contribute error to the operation simulated. 



operation performed, 
errors included, is 



u - H b zt d+Uf+Vc xy a?» cz ab xz- z™» x v b " , 



Specifically, the 
with all Pauli 
H a Z a lc ' - "a b zj b ^ ^ a ^ b 

where v x = 1 when i x = 4 or 5; w x — 1 when i x = 1, 3, 
or 5; u x = 1 when i x = 4 or 5 (with x being either a or 
6); and Uf = 1 when jf = 1, 2, or 5; u e = 1 when j e = 0, 
3, or 4; u c = 1 when i c = 2, or 3; Ud = 1 when id = 4, or 
5; v c = 1 when i c = 1, 3, or 5; and = 1 when id = 1, 
3, or 5. All the exponents are 0, otherwise. 



